Experimental estimation of one-parameter qubit gates in the presence of phase diffusion 
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We address estimation of one-parameter qubit gates in the presence of phase diffusion. We evaluate the 
ultimate quantum limits to precision, seek for optimal probes and measurements, and demonstrate an optimal 
estimation scheme for polarization encoded optical qubits. An adaptive method to achieve optimal estimation 
in any working regime is also analyzed in details and experimentally implemented. 
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I. INTRODUCTION 

Suppose you are given a black box which operates on 
qubits. The reconstruction of the corresponding quantum op- 
erations OIH Ht] is critical to verify its actions as a quantum 
logic gate J3l as well as to characterize decoherence processes 
J5D . Let us consider the case when the action of the device is 
described by the unitary = cxp{— icr n (j)} and corresponds 
to a phase-shift (rotation) <fr about a known axis n. This is 
the simplest operation on a qubit and realizes a one-parameter 
logical gate, which allows, combined with the Hadamard gate, 
the transformation of any qubit state into another. The quan- 
tum characterization U @, H H, [l(J [U d of this kind of 
devices is of interest for quantum information processing and 
amounts to estimate the phase-shift by measuring a suitable 
observable at the output. In the following, we fix the refer- 
ence frame and assume, without loss of generality, rotations 
U<j, = cxp{— io- z (j)} about the z axis. 

In ideal conditions the estimation of the phase-shift con- 
sists of preparing a qubit in a known pure state g = \ipo)(4>o\ 
and then performing a suitable measurements on the (pure) 
shifted state g^ = U^gUt, In a realistic implementation, 
however, the propagation of a qubit is unavoidably accompa- 
nied by some noise, which influences the estimation scheme 
and usually degrades the overall precision. In this paper we 
address estimation in the presence of the most detrimental 
kind of noise for a phase gate, i.e. non dissipative phase noise, 
which destroys the off diagonal elements of the density matrix 
and thus the information on the imposed phase-shift. We eval- 
uate the ultimate quantum limits to precision in the presence 
of noise and determine both the optimal preparation of the 
probe qubit and the optimal measurement to be performed at 
the output. The optimal estimation scheme is then experimen- 
tally demonstrated for polarization encoded optical qubits, to- 
gether with an adaptive method to achieve optimal estimation 
in any working regime. 



The paper is structured as follows. In section|II]we describe 
the system under investigation focusing on the ultimate limits 
to precision of phase-shift estimation in the presence of phase 
noise. The quantum Cramer-Rao limit as well as the optimal 
quantum estimator are explicitly given. SectionHIIl introduces 
a realistic scenario for qubit phase-shift estimation: an estima- 
tion scheme based on spin measurements is described in detail 
and the corresponding Fisher information is derived. In sec- 
tion|IV]we turn our attention to optical qubit systems and de- 
scribe phase-shift estimation under two different approaches: 
the inversion method and the Bayesian estimation. We also 
investigate numerically the robustness of the Bayesian analy- 
sis w.r.t. the inversion method in the non asymptotic regime, 
i.e., in the case of small amounts of data. SectionlVlis devoted 
to the experimental demonstration of the optimal scheme for 
phase-shift estimation, and of the adaptive method to achieve 
the quantum Cramer-Rao bound in any working regime. Sec- 
tion[VT]closes the paper with some concluding remarks. 

II. PHASE-SHIFT ESTIMATION IN QUBIT SYSTEMS 

The measurement scheme we are going to address is 
schematically depicted in Fig. Q] 
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FIG. 1 : Schematic diagram of the measurment scheme. 

One has a single qubit, initially prepared in the pure state 
g = |-0o)(V'o|, that undergoes an unknown phase-shift <f> im- 
posed by the unitary U^. Before being measured the shifted 
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state 00 = U^gU^ is degraded by a non dissipative phase 
noise occurring during the propagation. The effect of this kind 
of noise on the qubit density matrix, given a noise factor 7, can 
be described by the following Master equation (ME): 



£>0,a 2 = 7 L [ c^ + cr -]^ 



(1) 



where L[A]q^ A 2 = \ {[Aq^,^, ^} + [A Q^a 2 ^]}, and 
A 2 = jt/2 is, as we will see in the following lines, the ef- 
fective noise factor. Since L[er + (7_] and a z commute, we can 
focus the on the evolution of g, i.e., g A 2 = 7^[ cr + cr -]PA 2 - 
Upon writing g A -i in the eigenbasis of <j z , the ME leads 
to differential equations for the matrix elements g nm (t) = 
(n\g A 2\m), where g nm {t) = -57(11- m) 2 g nm (t) whose 
solutions read: 



,(o). 



(2) 



where g nm (0) are the initial density matrix elements. From 
Eq. d2|i it is clear that, whereas the diagonal elements are left 
unchanged by the evolution under ME (Q]), and, in turn, en- 
ergy is conserved, the off-diagonal ones are progressively de- 
stroyed. Finally, the solution of Eq. (0} is g^.A 2 = U^g A i U\. 
Since we can consider the noise factor A 2 as a fixed param- 
eter, in the following we will not write explicitely the depen- 
dence on it. 

It is worth noting that the same evolution as d2|i can be also 
obtained by the application of a random, zero-mean Gaussian- 
distributed phase-shift to a quantum state. Since the phase 
shift of an amount ip is described by the unitary operator U v = 
cxp(— i(p <j z ), we can write the state degraded by the Gaussian 
phase noise as follows: 
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e- lv{n - m) g nm (0)\n)(m\ (4) 
A ^ n - m >* Qnm (0)\n)(m\, (5) 



which is the same as in Eq. ©. This point will be useful for 
the experimental demonstration. 

The goal of an estimation problem is not only retrieve the 
actual value of the unknown parameter (the phase shift in our 
case), but obtain this information with the minimum uncer- 
tainty. The ultimate limit to the precision one can reach, i.e., 
the minimum variance, is given by the quantum Cramer-Rao 

bound d Hd d: 



Var[fl = [NH (<!>)]- 



(6) 



where N is the number of measurements and H(<p) is the 
quantum Fisher information (QFI). If we choose a pure probe 
state, g — \ip ) (ip 1 , then the QFI equals four times the fluctu- 
ation of <7 Z i.e one simply has H = 4(1 — (ipplcrMo) 2 ). More 
in general, for mixed states, H(<f>) can be written as IU7I1 : 



H{4>) = 2 Y, 



(A n — A m ) 2 

An + A m 



\{ip m {4>)\d^ n {<j>))^ 



(7) 



where the |V>n (</>)) 's are the eigenvectors of the state q$ = 
U(j,gU\ and A„ the correspoding eigenvalues. 

The spectral decomposition of g^ reads as follows 

n n 

\ip n ) being the eigenvectors of the initial state. If we decom- 
pose \ip n (4>)) in the standard basis as follows: 



\lpn(<f>)) = U^lpn) = E/^y"V w fc|fc), 



(9) 



then, by substituting (O into the eigenvalues equation 

g<p\tl>n{4>)} = A„|'0„(0)), after some algebra we obtain: 



A 2 (n-fc) 2 



qk 



\r qn Vn. 



(10) 



Moreover, since \d^tp n ((j))) 
[seeEq. ©]: 



*Sfc^ r ™fc e \k), we have 
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(ii) 



Finally, given A„ and r n k, we can evaluate the QFI, which is, 
as expected for a shift parameter, independent of <f>. 

Generally, the calculation of the eigenvalues A„ and the co- 
efficients r n k of the standard basis decomposition (0 is a dif- 
ficult task, which can be performed by means of numerical 
analysis. However, in the case of qubit systems, all the calcu- 
lations can be carried out analytically, as we are going to show 
in the following. 
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FIG. 2: Effect of the phase destroying ME on the Bloch sphere: the 
z component is left unchanged while x and y ones are scaled by a 

— A 2 

factor e . See the text for details. 



Upon writing the initial qubit state in the Bloch sphere rep- 
resentation 



1 + r ■ a 



(12) 
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where 1 is the 2x2 identity matrix, r = (r x , r y , r z ), \r\ 2 < 1, 
and a = {<J x ,Vy,Vz) is the Pauli matrices vector, then its 
evolution under the action of the ME (Q} can be reduced to the 
following transformation of the Bloch vector r lfl8ll : 



yxi r y> r z) 



I -A 2 -A 2 



0, 



(13) 



i.e., the Bloch sphere is deformed is such a way that the z 
component is left unchanged while x and y ones are scaled by 
a factor e~ A (see Fig. [2]). Now, due to symmetry consider- 
ations, without lack of generality we will focus our attention 
on the pure state with Bloch vector: 



r = (sin20,O,cos20), 



(14) 



29 being the azimuthal angle (9 = and 9 = tt/2 correspond 
to the north and south poles of the Bloch sphere, respectively). 

In the density matrix representation (choosing the a z eigen- 
vectors basis), after the phase-noise evolution, we have: 



cos 



e A cos 9 sin I 



The two eigenvalues are: 



1 ± 



1 + /(fl,A 2 ) 
V2 



(15) 



(16) 



with f(9, A 2 ) = ^e- 2A2 + (1 - e- 2A2 )cos46», and the cor- 
responding eigenvectors read: 



±{g±(e,A 2 ) 
z± \ i 

= ±[g ± (9,A 2 )\ + l/2) 



(tp+ltp-) = 0, and: 

g±(9,A 2 ) = cos20±/(0,A : 



1/2)], 



2 sin 26») 



Z± = + [g±(9, A 2 )] 2 



(17) 
(18) 



(19a) 
(19b) 



Substituting the previous equations into Eq. © we obtain (re- 
member that for qubit systems k = ±1/2): 



H(9,A Z ) = e- 2A sin^26», 



(20) 



which reaches the maximum for 9 = ir/4: the best states 
for phase estimation also in the presence of phase noise are 
the equatorial ones, i.e., the states laying in x-y plane of the 
Bloch sphere. 

Since the Bures metrics, and then, the Bures distance 
between states, is proportional to 
the QFI IU7[|26I1 . this result can be easily understood under a 
geometrical point of view. If we choose two states, one equa- 
torial and the other not, and shift them by the same amount 
(f>, then the distance on the Bloch sphere between the initial 
states and the shifted counterparts is larger for the equatorial 
states (see Fig.0 which, n turn, allow better estimation. 
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FIG. 3: The distance ab on the Bloch sphere between two equatorial 
states separated by a phase shift <f> is larger than the distance a'b' 
between non equatorial states separated by the same amount 4>. The 
same conclusion holds in the presence of phase noise: best estimation 
is achieved involving equatorial states. 



The optimal quantum estimator leading to < f2Qb can be writ- 
ten as lfl7h : 



H(9,A 2 )' 

where we introduced the symmetric logarithmic derivative: 

LrhPd> + o^La, 



(21) 



(22) 



and is the state of the qubit initially prepared in a pure 
state with Bloch vector r given by dT3l , shifted by the appli- 
cation of the unitary transformation = cxp(— ifyo+oJ) = 
cxp[— |0((T Z + I)] and degraded by the evolution through the 
noisy environment: 



e i0 A cog q gm 



e ^ A " cos 9 sin 
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One finds: 
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(23) 



(24) 



where X±, g± = g± (9, A 2 ) and Z± are given in ( fT6b and ( fT9l ), 
respectively. If we choose 9 = 7r/4, Eq. (|24T i reduces to: 



L = ie' A ' (cr+e** 



cr_e 



-id) 



(25) 



III. PHASE-SHIFT ESTIMATION BY SPIN 
MEASUREMENT 

Let us now consider a realistic scenario where, in order to 
estimate 4>, we measure the spin in a generic direction in the 
plane, i.e. the observable 

©a = &x cos a + a y sm a i (26) 

whose eigenvectors |S±(a)), 9 Q |S-|-(a)) = ±|E±(a)), are: 
|E ± (a)) = i=( e --| + l/2)±|-l/2)). (27) 
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For the qubit of Eq. d23T l. the probabilities to obtain the out- 
comes ±1 given the phase shift read: 



Pa*(±1\ 



Tr[|E±(a))(E±(a 
1 



lie 



—A* 



cos(a — 0) sin 29 



and the expectation value is: 

(6 a ) = Tr[6 Q q^)] = e~ A2 cos(a - 0) sin 29. 
The corresponding Fisher information turns out to be: 

P(0,A 2 ) = iM£#)[d lnP A2 (fc|0)] 2 
sin 2 (a- 0) sin 2 29 



fe=±i 



-2A' 



l- e - 2A2 cos 2 (a-0) sin 2 29' 



(28) 
(29) 

(30) 

(31) 
(32) 



which is plotted in Fig.|4]in the case of equatorial probe states 
(8 = it/ 4) as a function of 5 — a — and different values of 
A 2 . 
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FIG. 4: Plot of the Fisher information F{<f>, A ) in the case of equa- 
torial qubit probe states (9 = 7r/4) as a function of 8 = a — and 
different values of A 2 . 

As we can see, in the absence of noise (A 2 = 0) one has 
F = 1, Va, 0, i.e., the Fisher information is equal to the QFI 
H in Eq. d20b |27[ I2H1 . When noise affects the propagation, 
the maximum of F, that corresponds to the QFI H, is achieved 
for 5 = a. — <t> = tt/2, while goes to zero as a — = k tt, 
k € N. These results can be better understood considering 
the sensitivity of the measurement (actually this is the square 
of the sensitivity): 



S = 



Var[6 Q ] _ 1 - (0 Q ) 



(5 (e Q )) 2 (5 (e Q )) 2: 



(33) 



that is the ratio between the fluctuations of (O a ) and how 
(Go.) varies with respect to 0; ^/~S represents the smallest 
change of that can be detected with our measurement (up 
to the statistical scaling, of course). In the present case: 



S{<P, A 2 ) = 



1 - e 



-2A 2 



cos 2 (a - 0) sin 2 29 



e- 2A2 sin 2 (a-(/))sm 2 26l ' 



(34) 



which is just the inverse of Eq. (1321 ): the maximum of F (max- 
imum information) corresponds to the case of maximum sen- 
sitivity (minimum of S). If A 2 =0 one finds that Var[9 Q ] 



and (d^(Q a )) are always equal, no matter the values of a 
and <f>. When noise is acting, the maximum of F at S = tt/2 
corresponds to the minimum of 6> (</>, A 2 ), this fact can be also 
understood by geometrical means addressing the special case 
of a = (6o = <J X ). In this case the result of the measure- 
ment carried out out on the probe is just the projection onto 
the x-axis: for a fixed change d<p, the change of (0q) at = 
is smaller than the one at = tt/2 (see Pig . [5j ■ 

= tt/2 



degraded 
state 




FIG. 5: Sensitivity of a x measurement to a small change d(f> in the 
phase-shift. The projections onto x-axis are the expectations (Oo), 
with Bo = Ox- For a fixed change d(j>, the change of (Oo) at cj> — 
is smaller than the one at cj> = tt/2. For the sake of simplicity we 
sketch only the upper right quarter of the Bloch sphere x-y equatorial 
plane. 

The Fisher information depends on the actual, unknown 
value of 0. However, we can perform an adaptive, two-step 
method to achieve the QFI. During the first step, we use a 
small amount of data to obtain a rough estimate of the phase 
shift, then, at the second step, we tune 9 Q according to the 
transformation a — > + 7r/2: the (eventual) repetition of 
these two steps allows to reach the QFI limit. The same result 
can be obtained by fixing the measurement at a chosen a and 
tuning the probe state by applying a suitable rotation. 



IV. PHASE-SHIFT ESTIMATION FOR POLARIZATION 
ENCODED OPTICAL QUBITS 

In our proposal the qubit state corresponds to the polariza- 
tion degree of freedom of a coherent state. We refer to \H) and 
\V) respectively as horizontal and vertical polarization. Ini- 
tially, we set the polarization at |+) = -y=(\H) + \ V)) (equa- 
torial state) and apply the phase shift with eventual Gaussian 
noise (as we have shown in section|ll]this is equivalent to a non 
dissipative phase noise). Then we measure a x , which corre- 
sponds to inserting a half wave plate (HWP) at 22.5° in front 
of a polarizing beam splitter (PBS) and to record the number 
of counts n + and n_ at the two outputs, which correspond to 
the outcomes +1 and —1, respectively. Once a data sample 
X has been acquired, one needs to give an estimation of the 
phase shift. Here we consider two strategies: one based on 
the inversion of the probability functions, whose uncertainties 
come from the error propagation; the other applies the Bayes 
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theorem to retrieve a prior distribution of phase shift, knowing 
this distribution we can calculate all the desired moments. 



A. Estimation by inversion 

If we send M copies of our state, then we can write: 

n ± 



Var(<« 



P A2 (±1|0) 



n+ + ri- 



(35) 



and, by inverting Eq. (1291 ), after some algebra we obtain the 
following expression for the phase shift estimator: 



n+ - n- a : 
arccos [ e 



The uncertainty in the estimation is thus given by: 



(36) 



Var[0i, 



<7 2 (n+) + 



9n_ 



a\n-), (37) 



where a 2 (n± ) are the fluctuations of the numbers of outcomes 
n±. It is worth noting that at each shot, i.e., for each sent 
copy of the coherent state, the detected number of photons 
fluctuates according to the laws of quantum mechanics. 
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FIG. 6: Variances from simulated Monte Carlo experiments per- 
formed with M = 60 copies of the input state and A = 0.34. The 
plots on the top refer to the inversion method, the bottom ones to 
Bayesian analysis. For the plots on the left we chose a coherent state 
with an average number of photons equal to n = 2, for those on 
the right we set n — 12. The solid line is the Cramer-Rao bound 
l/(FMn). The dashed horizontal line is the quantum Cramer-Rao 
bound 1/ (H Mn). Notice the asymptotic optimality of the Bayesian 
estimation. 



C. Monte Carlo simulated experiments 



B. Bayesian estimation 



Given the data sample X = {x\, X2, ■ ■ ■ , xjy}, where Xk € 
{ — 1, +1}, Vfc, we can define the sample probability: 

N 

P(X\<f>) = Y[PA*(xk\<t>), (38) 

fc=i 

that is the probability to obtain the whole data sample X given 
the unknown phase <f>. By means of the Bayes theorem one can 
write the a posteriori probability l3lll : 

i N r 

P(0\X) = - l[P A 2(x k \4>), Af = /#P(^|X), (39) 
JV fc=i J * 

where $ is the parameter space. The probability ( |39l ) is the 
probability distribution of <f> given the data sample X. The 
Bayesian estimator are thus obtained as: 

<h= I d<t><t>P{<j>\X), (40) 
j<t> 

Var[0 B ] = [d<t>(cj>-<j> B fP{<t>\X). (41) 

Bayesian estimators are known to be asymptotically optimal, 
namely, they allow one to achieve the Cramer-Rao bound as 
the size of the data sample increases lfl7i l29i l30l l3~ll l32tl . 



Fig. [6] shows the variances of the estimated phases from 
Monte Carlo simulated experiments as function of the un- 
known phase shift. The simulations have been performed 
sending M = 60 copies of the input coherent state and choos- 
ing A = 0.34. The plots on the top refer to the inversion 
method, the bottom ones to Bayesian analysis. We inves- 
tigated two different cases: low energy (average number of 
photons equal to n = 2, left plots), and high energy (n = 12, 
rights plots). The solid line is the Cramer-Rao bound obtained 
from Eq. d34l ). while the dashed horizontal line is the quantum 
Cramer-Rao bound both rescaled by the effective number of 
measurements N = M n. We can see that in the low energy 
regime there is a more evident deviation from the expected be- 
havior (solid line) at the extremes of the phase interval (0, tt). 
This is the manifestation of a systematic error due to phase 
window sampling, i.e. the tails of the Gaussian distributions, 
simulating the phase-noise, is truncated at the boundary and 
refolded inside the interval. Notice the asymptotic optimal- 
ity of the Bayesian estimation when the energy, and, thus, 
the number of events increases. It is also useful to observe 
how the fluctuations of the variances obtained by using the 
Bayesian method are less than the other ones, which actually 
depend on the fluctuations of the average number of photons. 
A further investigation of Fig.|6]leads us to conclude that both 
the methods are robust over a large interval of phases: for the 
considered examples the variance is almost constant over the 
phase interval (0.5, 2.5). 
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V. EXPERIMENTAL SETUP AND RESULTS 

The experimental demonstration of our scheme is based on 
a KDP crystal, that allows both the manipulation of the opti- 
cal qubit polarization and the simulation of a phase diffusion 
environment. In Fig. [7] we sketch the experimental setup. 



HV 




FIG. 7: Sketch of the experimental setup to estimate the phase of 
a polarization qubit in the presence of non dissipative noise. The 
qubit state is set to |+) by inserting a half-wave plate (HWP1) along 
the path of the He:Ne laser output, while a KDP crystal, driven by 
a high voltage generator (HV), is used both to add the phase shift 
4> and to simulate the noise. A second half-wave plate @ 22.5° 
(HWP2) together with a polarizing beam splitter (PBS) implement 
the measurement (the couplers, C, and the multi-mode fibers, MF, 
are used to bring the signals to the avalanche photodiodes, APD). Fi- 
nally, the contemporary counts Ni g and N2 g coming from the APD1 
and APD2 and from the APDg, i.e., the gate channel signal com- 
ing from the beam splitter (BS) placed in front the laser source, are 
acquired by a PC module. 

A linearly polarized He:Ne laser (Thorlabs HRP120) gen- 
erates coherent states that are prepared in the initial qubit state 
|+) = -j=(\H) + \V)) by means of a first half-waveplate 
(HWP1). Then the qubits pass through a KDP crystal driven 
by a stabilized high voltage generator (HV), up to a maximum 
of 6 kV. The application of an electrostatic field E z along the 
z axis of the KDP (that is also the direction of propagation 
of the qubit), a phase shift <j>{E z ) is introduced between the 
qubit polarization components along x and y. The crystal is 
oriented with the x axis parallel to the horizontal polarization 
and tilted around the vertical axis y in order to avoid the mul- 
tiple internal reflections superposed with the main beam, that 
are removed by a pin hole. 

The detection system consists of a second half-waveplate 
(HWP2), a polarizing beam splitter (PBS), absorption filters 
and three detectors. The HWP2 plate is set at 22.5° with re- 
spect to the horizontal axis and the PBS selects the polariza- 
tion. Light signals are focused into multi-mode fibers (MF) 
and sent to the detectors APD1 and APD2. The detectors are 
single photon counting modules (SPCM) based on avalanche 
photodiodes operated in Geiger mode with passive quenching. 
For the coincidence counts an electronic circuit based on AND 
gates is used. A 50:50 Beam Splitter (BS) is placed before the 
HWP1 plate to add a gate channel (g) for the counting mea- 
surement. The aim is to make acquisition with a low number 



of counts maintaining an high signal to noise ratio Ni/Ni dc 
(i = 1, 2). This is the reason for the coincidence counting, 
according to the formula (i = 1, 2): 

N ig = {Ni + N udc )N g At = N t N g At + N i>dc N g At 

= Nigjrue + N ig ^dc, (42) 

where Ni g , Ni and N g are respectively the coincidence and di- 
rect counts on the detectors APDi and APDg, Ni g ^ c and Ni^c 
are respectively the coincidence and direct dark counts, and 
At is the coincidence time window of the electronic count- 
ing module (At = 90 ns). Ni g could be as low as we want 
by selecting the time window of the acquisition and the ratio 
Nig^true/Nig^c = Ni/Ni^ c remain constant. We measure 
about 10 5 counts/s on the detector APDg and a maximum of 
about 9 x 10 4 counts/s on the other two but we can adjust 
these direct counts changing filters or the detector voltages 
above breakdown. The dark counts are below 200 counts/s. 
Each acquisition is taken with a 10 ms time window in order 
to have n ~ 10 in the contemporary counts. Note that the sig- 
nal to noise ratio remain ~ 10 3 . The rising and falling time 
of the high voltage generator are measured around 200 ms so 
the waiting time between two subsequent measurement is set 
at 240 ms. 

We made M=60 acquisitions for each phase-shift <f>. Notice 
that each acquisition corresponds to a different random phase 
distributed according to a Gaussian distribution centered on <fi 
and with standard deviation cr 2 (</>) = A 2 /2 [see Eq. (|3)]. In 
our implementation for each acquisition we send a different 
voltage to the KDP crystal according to a proper calibration 
curve <p(V). 

The results are shown in Fig. [8] where we plot the theo- 
retical Cramer-Rao bound (solid line) and the experimental 
variances obtained using the inversion method (red squares) 
and the Bayesian analysis ( blue circles). The Cramer-Rao 
bound is quickly reached by Bayesian estimation, even for a 
relatively small number of measurements (N = Mn in the 
figure, with M = 60 and n ~ 10). 



A. Adaptive method for optimality 

As we have seen Bayesian estimation allows us to reach the 
Cramer-Rao bound for any orientation of the measurement an- 
gle. On the other hand, as we pointed out in the section Hill 
optimal estimation, i.e. with precision at quantum Cramer- 
Rao limit is achieved only if <fi = ir/2. In this Section we 
show that ultimate precision can be always achieved upon the 
application of an adaptive method l32l mi l34tl . In the present 
case, we proceed as follows: i) we start our estimation obtain- 
ing a first value fa of the unknown phase; ii) we set a of the 
measured observable (l26t to the value cti — ir/2 + fa (this 
can be also obtained by adding the phase shift a[ = it/ 2 — fa 
to the qubit itself) so the new phase is near to 7r/2. One can 
also repeat many times this procedure in order to refine the 
estimation. However, as one can see in Figs.[6]and[8] the op- 
timal condition is almost achieved after the first step and the 
variance is nearly constant around the optimality region. We 
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FIG. 8: (Color online) Variances from experiments performed with 
M — 60 copies of the input state and, from top to bottom, A = 
0.13 ±0.02, A = 0.24 ±0.03 and A = 0.48 ±0.05. (Red) Squares 
refer to inversion method, (blue) circles to Bayesian analysis. We 
also have, from top to bottom, n = 9.83 ± 0.93, n = 11.06 ± 
0.43 and n = 9.78 ± 0.95. The solid line is the Cramer-Rao bound 
l/[FMn). 



qubit in the presence of phase diffusion. Our scheme is based 
on polarization measurement assisted by Bayesian estimation 
and allows to achieve the ultimate quantum bound to precision 
using a limited number of measurements. In turn, Bayes es- 
timator is known to be asymptotically unbiased, but for prac- 
tical implementation is of interest to evaluate quantitatively 
how many measurements are needed to achieve the asymp- 
totic region. Our results indicate that Bayesian inference rep- 
resents a useful tool for phase estimation. An adaptive method 
to achieve optimal estimation in any working regime, i.e. for 
any value of the unknown phase-shift, has been also analyzed 
in details and experimentally demonstrated. As a future per- 
spective we foresee the possibility to investigate the effects 
of different kind of noises and to employ entangled probes to 
increases the overall stability of the estimation procedure, a 
results which has been theoretically established for the ideal 
(i.e without noise) case l3lll . 



step Var(0) 






1 4.21 x 10" 


3 


0.33 ± 0.06 


2 2.49 x 10" 


-3 


0.22 ± 0.05 


3 2.52 x 10" 


-3 


0.19 ±0.05 


4 2.59 x 10" 


-3 


0.20 ± 0.05 



TABLE I: Phase estimation using the adaptive method (see the text 
for details). The actual value of the phase is cj> — 0.17 ± 0.01. The 
amplitude of the phase noise is A = 0.46±0.06 and the energy of the 
coherent state is n — 10.97 ± 0.67. For each step we used M — 55 
repetitions for the estimation. For these values of the parameters the 
quantum Cramer-Rao bound, evaluated from Eq. J20b and rescaled 
by the number of measurments N = Mn, is equal to Var(</>) ~ 
2.52 x 10" 3 



report the results in Table[T] where one can read the step num- 
ber of the adaptive method, the corresponding variance and 
the estimated phase. Note that already at the second step op- 
timal estimation has been achieved, that is the experimental 
variances are close to the optimal variance given by the quan- 
tum Cramer-Rao bound. 

VI. CONCLUSIONS 

We have demonstrated an optimal estimation scheme for 
the phase-shift imposed to a polarization enconded optical 
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